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THEtRY  cf  shallow  donor  states  with  a strong  central  cell  perturbation 


Joan  R.  Cnffrcy 


Dept,  of  Physics,  University  of  California,  San  Diegc 
La  Jolla,  California  92093 


A new  theory  for  shallow  donor  states  with  a strong  central  cell 
perturbation  is  presented.  The  impurity  wave  function  is  expanded  in 
terms  of  generalized  Wannier  functions,  equivalent  to  the  eigenfunctions 
of  the  periodic  lattice  with  the  strong  central  cell  perturbing  potential, 
"he  single  unknown  parameter  of  this  theory  can  be  fitted  by  comparison 
with  the  observed  binding  energy.  The  theory  is  applied  to  the  impurity 
states  of  As,  p and  Sb  in  Silicon.  The  remi  contact  constants  for  the 
Si^G  lattice  nuclei  surrounding  the  impurity  are  calculated  and  compared 
te  the  E1ID0R  experiments  of  Hale  and  Pieher  and  to  recent  calculations, 
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THEORY  CF  SHALLOW  TNCR  STATES  WITH  A STRONG  CENTRAL  CELL  PERTURBATION 
I.  INTRODUCTION 

Kuch  theoretical  work  has  been  done  on  the  calculation  of  impurity 
energy  levels  and  wave  functions.  The  shallow  donor  energy  levels  in 
Silicon  were  qualitatively  understood  very  early  in  terms  of  a hydrogenic 
model  called  the  Effective  Mass  Theory  (EMT)"*,  in  which  the  impurity  wave 
function  is  expanded  in  terms  of  the  Bloch  waves  of  the  perfect  crystal. 

For  sufficiently  localized  perturbing  potentials,  the  Koster-Slater 

O' 

theory  , in  which  the  impurity  wave  function  is  expanded  in  terms  of 
.’annier  functions-5,  is  often  used.  lecently  pseudopotential  theory 
has  been  adopted  to  calculate  the  impurity  energy  levels  and  wave  functions. 
Pantelides  and  Sah  , as  wel 1 as  Schechter?-0,  have  reformulated  the  FIT 
in  terms  of  a pseudopotential  scheme.  However,  in  both  cases,  the  use 
of  the  E?T  was  logically  not  justified  because,  although  their  pseudo- 
potentials were  weak,  they  were  not  slowly  varying  on  the  scale  of  the 
lattice  cell  dimension. 

Appro yimately  simultaneously  with  this  research,  a detailed, 

r.any-band  calculation  of  impurity  states  in  Silicon  was  undertaken  by 

Ivey  and  Vieher  (i-K)  Their  calculation  was  very  thorough  and  was 

the  most  comprehensive  analysis  of  the  wave  functions  for  shallow  dcnor 

states  in  Silicon  to  date.  Their  results  were  in  much  better  agreement 

with  the  experimental  data,  from  the  ENDO.R  experiments  of  Hale  and 
11  12 

Fieher  , than  those  of  any  previous  work.  I-K  compared  their  results 

with  those  given  by  the  FIT  and  concluded  that  the  calculation  of  these 
shallc”  dcnor  i-purity  "ave  functions  is  outside  the  realm  of  the  E!T. 

Since  we  have  just  completed  a calculation  of  these  donor  wave  functions, 


-2- 


using  a generalized  fern  of  the  Ef!T,  we  will  conclude  with  sene  comments 
on  their  method  and  conclusions. 

We  believe  that  the  present  approach  succeeds  for  the  first  time 
in  yielding  a logically  correct  one-band  theory  of  shallow  impurity 
states  with  a strong  central  cell  perturbation.  The  effects  of  both 
the  deviations  of  the  potential  from  Coulorbic  form  and  the  important 
interband  effects  produced  by  the  central  cell  potential  are  describable 
by  generalized  Wannier  functions  (GWF)  belonging  to  a single  band.^ 

From  a theoretical  point  the  shallow  impurity  states  in  Silicon 
are  not  a particularly  ideal  case  for  the  application  of  this  theory 
since,  due  v the  small  energy  gap  between  the  valence  and  conduction 
bands,  the  Wannier  functions  are  not  very  highly  localized.  The  reason 
why  these  systems  were  nevertheless  chosen  was  that  all  the  necessary/ 
information  about  the  structure  of  the  conduction  bands  and  the  impurity 
states  was  -mown.  Another  class  of  systems  that  might  be  successfully 
treated  by  this  approach  are  color  centers  in  alkalai  hal^ides. 


II. 


GENT  lr..  ’"Th  MEl  FUNCTION  THEORY  C?  IMPURITY 

applications  tc  donor  stater  in  sms  on 

A.  Formulation  of  the  Theorv 


STATES  WITH 


'•'.’hen  a shallow  donor  impurity  is  introduced  into  a perfect  lattice, 
the  Hamiltonian  for  the  electrons  may  be  written  as 


v = u 


v , 


(D 


where  H0  is  the  Hamiltonian  of  the  perfect  lattice  and  U is  the  perturbation 
due  to  the  impurity,  which  has  the  asymptotic  form. 


u(?) 


Kr 


(2) 
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( K is  the  static  dielectric  constant  of  the  host  lattice.)  The  wave 
functions  satisfy  the  differential  equation, 


(H  - )W  = C. 


(3) 


The  ground  state  wave  function  ray  be  expanded  in  terms  of  the 
perfect  lattice  Vannier  functions,  a°^Q,  as  in  the  Koster-Slater 
theory^"  - f 


IV 


I c 


n,5  n,s 


•here  n is  s band  index  3nd  s is  a site  index, 
the  set  cf  difference  equations, 

? ^ ^ an', s^  cn', s'  _ * cn.et  * 


(U) 


"he  coefficients  satisf-, 


(5) 


If  the  potential  is  a smoothly  varying  function  cf  r,  the  matrix  elements 
cf  H will  be  essentially  mere  unless  the  band  indices  are  the  same,  in 
- hich  case  (S’)  reduces  to  a one-band  set  cf  equations.  However,  in 
the  shsl'r"  donor  impurity  problem,  the  potential  is  net  smoothly  varyin^ 
r.r.  the  central  region  sc  that  solution  of  the  set  cf  equation*-  {■■) 
becomes  a very  difficult  many-band  problem. 

In  order  to  circumvent  this  difficultv  we  expand  the  wave  function  ir 
terms  of  an  appropriate  set  of  generalised  'larmier  functions  (G'.vF),  which 
are  constructed  so  33  t(  insure  that  the  matrix  elements  of  H between 
these  functions  is  negligible  unless  their  band  indices  are  the  same.  V.’e 
write  the  total  potential  as  the  sum  of  two  terms 

U(r)  = U,(?)  + U2(r) 


J 


(6) 


’•here  T^"?)  is  a smoothly  varying  function  of  "r  for  all  r,  which  reduces 
to  the  fern  ;iven  by  (2)  in  the  asymptotic  region  and  is  smoothly 
continued  throughout  the  central  region  in  an  arbitrary  ray.  3y 
smoothly  varying  T_e  mean  that  UgC?)  is  essentially  constant  over  the 
range  cf  a h'annier  function,  or  or  cur  requirements  on  U 2 ( r) , U 1 Or)  is 
necessarily  short  range. 

In  a manner  similar  to  that  described  in  a recent  paper  by  the 

..  1'.  ... 

a\’  ,hor  arc  ..  I.cnn  •,  re  can  construct;  one  ,,  "nxch  are 

equivalent  tc  the  eigenfunctions  of 


in  which  the  short  range  potential  U-)  is  included.  These  GMF  are 


orthonormal , 


(xivi  ~ Ng  cV  " Sn  v-/  c7  • 


They  have  the  same  degree  of  ? realization  as  the  perfect  lattice 
ar.nier  fun  stir  ns  an'3  a^orcoch  them  exponential  Iv  for  lattice  sites 
distant  f-sm  the  impurity  site.  Their  construction  insures  that  the 
matrix  elements  of  H-j  between  them  is  zero  unless  their  band  indices 
are  the  sa-e.  Since  Ug  is  required  to  be  smooth,  all  cf  its  off- 
diagonal  -atrix  elements  in  both  n and  s will  be  negligible. 

V.'e  nor  -■.•rite  the  ground  state  impurity  wave  function  as  a linear 
combination  of  the  GV.T  associated  with  the  lowest  band  (and  drop  the 


band  index) , 


W?)  ■ I c.  as(r)  , 


■•here  the  cs's  satisfy  the  set  cf  difference  equations 
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(“c  I -lasV  c„,  - 

5 *"*  **3 


do) 


In  the  region  far  from  the  irpurity  the  GY?  arc  the  sane  as  the  perfect 
lattice  Yannier  functions;  U 1 ( F)  is  zero;  and  l:2(r)  is  given  by  (2). 
Therefore,  the  Hamiltonian  ret six  element  con  ne  written  as 


(a,nias.)  = (a?lHcla°,)  - ll  8s>n»  - A.. 


(11) 


Kr, 


’’here  rc  is  the  location  of  the  nidpcint  of  the  sth  Yannier  function  and 
the  matrix  elerer.tr  A.  describing  the  effect  of  the  short  ranre  no  ten- 
y ore*,  ’V  -h?  r;  one  rr;,  ainkriov?n  •;  ’or.e . .c  a : -ire 


(12) 


^ r-  I “ ) 

' J-1 

eroert  f*r  Yannier  function?  near  the  i-’-urity . 

“ • A~  -tcation  uo  honor  States  in  Silicon 

• -.-ill  nov  consider  a shallcr  aonor  inpurity  ir;  Silicon.  The  four 


” ~v 


lien  band*  in  this  crystal  are  connected.  As  was  rho’-n  in  a oaner 


oy  hr  hr  , in  order  to  obtain  a localised  ret  of  antibenSin"  V.  annier 

> ~t  the  eigenfunction-  as  so  elated  vitb  there 

hands.  Inch  ’.  annier  function  i 1- belled  by  a site  index  i an-,  a bend 
i-.’.e.v  r,  "here  r.  refer-  to  one  cf  the  four  bonds  associated  with  the  a ter 
at  site 

The  ground  state  wave  function  transferee  according  to  the  smretric 
re-re-entation  A1  of  the  crystal  point  group  T.- . Hence,  those  expansion 
ccef**4.  cient?  "hid  are  e-uivrler.t  under  the  crer-tirn-  of  the  -roue  T,  are 
eniE'  . This  leads  one  to  define  a new  set  cf  erthonema?  basis  functions 
ar,  tier  are  proportional  tc  the  sur,  of  all  .’annier  functions  which  are 


ki 


, '=  6S,1 


(14) 


TY.p  mtri'  elements  of  Hq  between  the  perfect  lattice  h'annier  functions 
are  related  tc  the  Bloch  energies  E°  tj  and  were  determined  by  a fit  to 
the  Silicon  energy  conduction  bands. 

The  single  parameter  A was  deteniined  by  fitting  the  experimental 
binding  energy,  Eexp.  If  we  define 

Ks,s'  = (^lH0laO,)  - £_  8SfS,  - Sexp  5 , (15) 

Krs 

equation  (11)  ray  be  written  as 

~ ^^S,S'^CS'=  (16) 

where  A ray  oe  regarded  as  an  unknown  eigenvalue,  given  by  the  stationary 


>1 


I 

I 


' 

1 


! 
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expression 


A - 


I 

s,s' 


I 

s,s 


t 


GS  Is,s' CS' 


(17) 


which  is  a minimum  for  the  correct  cs's.  Because  of  computer  limitations 

we  were  only  able  to  determine  Li?  unknown  coefficients.  We  assumed  that 

the  rest  were  given  by  their  EKT  values.  For  the  impurities  Sb,  P,  and 

As,  the  parameters  A were  determined  to  be  .187,  .219,  and  .281  eV 

resnectively.  We  found  that  the  coefficients  for  the  Wannier  functions 

far  from  the  impurities  differed  very  little  from  their  EITT  values. 

29 

The  Ferri  contact  constants,  a.,  at  the  Si  nuclei  surrounding 
the  impurities  (?,  As,  Sb)  have  been  determined  from  R’JDOR  experiments 
performed  originally  by  G.  reher  1 and  most  recently  and  completely  by 
Hale  and  I'ieher11.  They  are  proportional  to  the  wave  function  density 
at  the  site11, 

IUI(R.)|2  = 1.523  x 10"2  I1--,  a. | X-3  , ( IS) 

J v 

where  a • is  expressed  in  HHz . When  the  crystal  orientation  in  the 
J 

magnetic  field  is  varied,  the  angular  dependence  of  the  ENDOR  signal 
allows  one  to  recognize  which  signals  belong  to  nuclei  on  sites  of 
the  types  <aaa>,  ^OOa'?,  and  <aab\  A knowledge  of  the  wave  function 
at  the  lattice  sites  enables  one  to  match  observed  experimental  lines 
to  specific  lattice  sites. 

In  the  calculation  of  the  Fermi  contact  constants  we  used  orthego- 
nalized,  directed  tetrahedral  orbitals  for  the  Wannier  functions.  Our 


j 


i 
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results,  along  with  those  of  I-K  and  the  experimental  values11,  are 

presented  in  Tables  I and  H . We  agreed  with  I-K  on  the  matching  of 

slightly  more  than  half  of  the  observed  lines  with  specific  lattice 
1 ft 

sites  . Agreement  with  I-K  for  sites  of  the  type  <aaa>  and  <00a^  was 
perfect.  For  sites  of  the  ether  type,  the  matching  procedure  is  more 
ambiguous,  and  we  have  less  confidence  in  the  assignments  offered  by 
either  I-K  or  by  us.  The  overall  agreement  of  our  values  with  experi- 
ment is  slightly  worse  than  that  of  I-K,  presumably  because  of  the 
relatively  poor  localization  of  the  Silicon  Wannier  functions. 


> 


I 

| 


in.  sc:e  comments  on  the  ivey-iheheh  method  and  conclusions 

Ive”  ar.d  Kieher1^  use  a one-step  iteration  to  calculate  the 
pseudo  impurity  wave  function,  which  satisfies  the  equation 

(Ho  + Up  - E)Wp  = 0,  (19) 

where  Up  is  the  pseudopotential  and  is  expanded  in  term’s  of  the 
pseudo  Bloch  waves  with  coefficients  An(k) . HQ  is  the  Silicon  lattice 
Hamiltonian,  which  is  highly  anisotropic.  To  lowest  approximation 
they  use  a function  which  has  spherical  symmetry.  To  allow 

for  the  anisotropy  of  the  wave  function  due  to  the  anisotropy  of  HQ 
and  to  obtain  the  components  of  the  wave  function  associated  with 
other  branches  of  the  conduction  and  valence  bands,  thev  then  solve 


(H 


o 


= - Up 


v 20) 


This  is  not  a consistent 
and  therefore  we  decided 
ar.  estimate  of  the  errors 


iteration  technique  for  obtaining  the  anisotropy, 


tc  do  a simple  calculation  in  order  to  obtain 
inherent  in  this  procedure.  Cn  the  basis  cf 


this  calculation  we  estimate  that  the  anisotropy  obtained  from  (20)  is 
up  tc  30“?  too  snail. 


e of  a Silicon  lattice  effective  Hamiltonian 


■opv 


•rave  function  satisfie 
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Thus  we  see  that  use  of  the  procedure  (20)  introduces  significant 
errors  into  those  parts  of  the  wave  function  T-'hich  are  associated  with  the 
two  lowest  conduction  bands,  fn  the  other  hand  one  can  verify  that  the 
I-”  equation  (20)  does  give  a gccd  appro  yin  r-tion  for  the  parts  of  UJ. 
associated  with  the  other  branches  of  the  conduction  and  valence  bands. 

In  their  abstract  I-!'  state  that  the  EH  restricts  itself  to  a 
single  band  and  to  conduction-band-minina  Bloch  functions.  We  would 
like  to  point  out  that:  (l)  When  an  energy  band  ninirrur.  is  near  a 

band  edge,  the  V!*T  rust  employ  an  extended  zone  schere  which,  in  a 
reduced  zone  schere,  le'ds  to  inclirlon  of  ~ov*e  than  one  b-nd . (2)  ?he 
1:1  assures  that  the  wave  functions  are  nr.de  ’r  of  -ave  ••'ackets  of 
iv.  ch  functions,  localised  .near  the  conduction  band  ninira.  This 
is,  in  fact,  confirmed  b;*  the  results  of  I-’’,  as  the--  find  that  the 
•'.■a  !r  e:-:c*  ' s ’ in  the  (er-tendel  rone/’  '1  contribute  t'-ei  ml” 

about  11"  to  the  Terri  contact  constants. 


There  are  two  level':  of  annre  xinaticn  in  the 


3 ~ . i . on  the 


-are  acc-’r-te,  level  of  anprc'i-roticn  the  leurier  ir-  nrferr,  1.(r),  of 

the  loch  "ave  coefficients  Apr)  (taken  with  respect  tc  the  conduction 
. 

o?no  nzninur,  kj)  satisfies  an  anisotropic  hydro  genic  t;/pe  differential 

c 

eguaticn.  "hr>~  Apr)  is  rea'  and  even  in  the  components  of  (T:-k^)  . Cn 
the  second,  less  accurate,  lever  cf  arr.ro xirstion  the  irnurlty  wave 
.uncticr  at  the  lattice  s-tes  ' is  •-ritten  ?s 

"jr  I ?,(V^  (g),  («) 
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l>-  ■ 


AUO  near  In  fact,  for  "115.000  this  is  net  a good  approximation  for 

J 


the  lattice  sites  rjhiflh  do  not  have  inversion  s'rr.etry,  'er.ee  tbit 
second  level  of  ;; cproxiration  cannot  be  expected  tc  hold  for  these 
sites . 


I-M  calculate  a function  F.  defined  by 

J 


(*£.)  = jf'  A(ic) , 


(27) 


vhieh  they  carps  re  T'ith  the  FT  function  F.(R-).  (The  superscript  j on 
the  -urn — tion  indicates  that  the  surnation  is  perf'-rred  ever  the  one- 
sixth  of  k-space  associated  with  the  j^h  conduction  band  minimum  and 
extendin'  through  the  first  tvo  ’Jriliouin  nones.)  The  - find  that,  for 
nan/  of  t::e  .attice  aites  .(h.'.ch  dc  not  }."ve  inversion  a ’r’.netr/,  . 

O 

la.,  a "ary.,  ir  aginar"  ~ri.  They  conclude  that  their  results  are  net 

carp  at  able  adth  the  FT.  In  fact,  the  main  inconsistency  is  with  the 
second,  les~  accurate,  leve"  of  the  fT  since  fer  rest  o?  t,he  lattice 

.■a 

ri  ter  vith  inversion  — retr-*  the  imaginary  rt,  of  is  sr.-ll.  The 
-••’in  exception  to  this  statement  is  IF  for  the  site  (0,0, c),  -hich  1-1 
report  to  have  an  imaginary  part  'daich  is  about  IQf  of  its  real  part. 
Since  we  were  surprised  by  this,  we  also  calculated  fer  this  site 
using  our  sve  function.  We  found  that  the  imaginary  part  was  about  20,1 
cf  its  ''er-  "art  rn  ' that  ro~t  of  it  (about  Hi')  --as  are  te  a real  ter~ 
in  A(k)  which  va-  odd  in  (S-k*),.,  and  which  arose  fror  the  odd  nsrt  cf 
the  pc central  matrix  e erer.t  ( ^ | |l|J) . Thus  this  'arge  imaginary  mart 
does  not  signify  a breakdown  of  the  FT. 

We  believe  that  a corps  risen  of  our  work  and  that  of  I-'',  is  useful 





■J 
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for  the  understanding  of  shallow  donor  states  and  specifics  }!*'■  to  the 
application  of  the  to  this  problem,  ’\Te  found  that  by  u^inr 
generalised  Cannier  functions  and  by  matching  a good  inner  region 
solution  to  the  asymptotic  TI TP  solution  (which  only  u-e  ' the  f j rs t 
level  of  appro miration) , vTe  ’.rere  able  to  calculate  quantitative  impurity 
state  wave  functions,  including  strong  central  cell  effects. 
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